Stripe formation in electron-doped cuprates 
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We investigate the formation of charge domain walls in an electron-doped extended Hubbard 
model for the superconducting cuprates. Within an unrestricted Hartree-Fock approach, extended 
by the introduction of slave-bosons to obtain a more proper treatment of strong correlations, we 
demonstrate the occurrence of stripes in the (1,1) and (1,-1) directions having one doped electron 
per stripe site. The different filling, direction and width of these electron-doped stripes with respect 
to those obtained in the hole-doped systems have interesting observable consequences, which are 
discussed. 
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The issue of charge inhomogeneity in the cuprates and 
particularly in the underdoped region of their phase dia- 
gram is presently attracting a great attention both, from 
the theoretical and the experimental point of viewd. 

Some theoretical proposals are based on the occur- 
rence of stripe phases with the charge tending to order 
in incommensurate charge-density waves or in fluctuat- 
ing domain walls to explain both the awwnalous normal 
state properties and. superconductivityBa. In particu- 
lar it was proposed] that the stripe phase in the un- 
derdoped cuprates is ruled by a charge-driven critical 
line ending in a hidden quantum critical point (QCP) 
located around optimal doping. Although this QCP may 
be masked by (strong) pairing .due to the nearly singular 
attractive charge fluctuationscm the physics of the under- 
doped cuprates is strongly influenced by the tendency to 
spatial charge order. 

On the experimental side, the critical line for stripe 
formation has likely been detected in NQR0 and in mag- 
netotransport experimentajfl, where the suppression of 
superconductivity by strong magnetic fields both in hole- 
and electron-doped materials reveal a metal-to-insulator 
crossover line ending at T = near optimal doping. In 
the hole-doped compounds the insulating character of 
the low-temperature and low-doping phase is stronger at 
the specific doping 1/8, where commensurability effects 
between the charge ordering and the underlying lattice 
make the charge order particularly robustQ. This latter 
effect is a signature of spatial order and provides fur- 
ther evidence of the enhanced strenght of stripe ordering 
in various compounds at specific fillings (this is the so- 
called 1/8 phenomenon). This commensurability effect 
seems instead much less evident in the electron-doped 
materials. 

There is also an increasing and more direct experimen- 
tal evidence for the occurrence of (static or dynamic^ 
charge ordering in under- and optimally dopecLeupratesEJ, 
ranging from |EXAFSEj to neutron scattering^, and X- 
ray scatteringE3£3. 

The stripe instability was predicted theoretically in 
Rcf. [lJI within Hartree-Fock (HF) theory applied to the 
extended Hubbard modci-and confirmed by a number of 
subsequent investigational^. Only recently, however, the 



occurrence in the single-band Hubbard model of half- 
filled stripes (i.e. with half a hole per stripe site) in 
the (1,0) [or equivalently in the (0,1)] direction [Half- 
Filled- Vertical stripes (HFVS)] was demonstrated within 
a generalized HF approach including slave bosons to 
suitably handle the strong local electron-electron inter- 
action. This achievement remarkably showed that (a 
proper treatment of) strong correlations and long-range 
interactions are crucial in determining a specific stripe 
texture similar to the one observed in the cuprates. 

The aim of the present work is to apply the same 
joynt HF-slave-boson technique to theoretically investi- 
gate the stripe formation in a more realistic extended 
three-band Hubbard model in the case of electronic dop- 
ing. It is to be borne in mind that in any HF approach 
the solutions are self-consistently stabilized by locating 
the chemical potential in the middle of a gap. There- 
fore this approach is only suitable to describe insulating 
phases where the filling and/or the breaking of transla- 
tional symmetry gives rise to such gaps. In our case we 
therefore intend to describe by our method the 1 /8 com- 
mensurate systems where both the charge ordering and 
the (nearly) insulating character are substantial. 

— The model and the main outcomes. The starting 
point is the extended three-band Hubbard model 

H = E nf{di Pja +h.c.)+ ]r t^ipip^+h.c.) 

<ij>a <jj'><? 

+ A XX'''- + Ud Yl d h d ^ d li d n 
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where the first and second terms describe the nearest 
neighbor (nn) copper-oxygen and oxygen-oxygen hopping 
of holes in the d copper orbitals and in the p oxygen or- 
bitals respectively. A is the energy difference between 
p and d orbitals (the energy of the d orbitals is set to 
zero), while Ud {U p ) is the local Hubbard repulsion on 
the copper (oxygen) orbitals. The last term is the long- 
range interaction, with Vtj = Ve~ rij t Xs / (2ry), with ry 
being the distance between two sites, As is the screening 
length. V is set to V = Upde 1 ^ 3 in order to match the 
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Upd mi Cu-0 repulsion. In this term, the c operators 
represent d or p hole operators depending on the nature 
of the site. 

Following the procedure of Ref . |l^ we treat the strong 
local repulsion Ud by means of a slave-boson approach. 
The weaker strength of the repulsion on p orbitals and 
the smaller average hole density on oxygen sytes makes 
unnecessary to apply the same technique to the oxygen 
orbitals. Therefore the other interaction terms U p and V 
are treated by a standard HF decoupling. 

We summarize here the main results valid both for 
hole and electron doping, i) Like in the single-band 
Hubbard model, also in the three-band Hubbard model 
the major effect-pf the slave-boson treatment is to stabi- 
lize the stripesEJ with respect to isolated spin polarons. 
This stability mostly arises because, according to gen- 
eral argumentscj, the charge mobility tends to favor 
stripes with respect to isolated spin polarons. Then, since 
the Gutzwiller-like procedure implemented via the slave 
bosons allows the charge to delocalize without paying a 
much too large repulsive energy stripes are not artificially 
disfavored within our approach. The same is not true in 
standard HF treatments, where the double occupation of 
the sites is avoided via a strong antiferromagnetic polar- 
ization of the spins, which greatly suppresses the inter- 
site hopping, ii) For moderate to large values of Ud, the 
natural energy scale for charge excitations no longer is 
Ud, but rather the charge-transfer energy A becomes the 
principal parameter determining the stability of stripes 
with respect to an ordered lattice (Wigner crystal) of 
maximally separated spin polarons. The stripes are only 
stable below a critical value A c , depending on the other 
parameters of the model, which in any reasonable realis- 
tic case is much larger than the critical values obtained in 
the standard (i.e. without slave-bosons) HF technique. 

iii) For a quite broad range of realistic values of the 
parameters the HFVS and the Filled Diagonal stripes 
[i.e. with one hole per stripe site and running along the 
(1,1) or (1,-1) directions (FDS)] are the most stable or 
the only stable solutions. In particular, for hole dop- 
ing x — 1/8, while the standard HF treatment favors the 
FDS, it was remarkably found that the HFVS become the 
ground state when strong correlations are properly treated 
by slave bosons. The absolute stability of HFVS occurs 
even in the absence of long-range interactions. In this re- 
spect, the hole-doped three-band Hubbard model seems 
more favorable than the single-band one in giving rise to 
half- filled stripe textures of the type observed in cuprates. 

iv) The local repulsion on oxygen only plays a minor 
role, while v) the oxygen-oxygen hopping t pp enhances 
the charge mobility and therefore stabilizes the stripes 
with respect to the polaron Wigner crystal (POLWC). 
As far as the effects of the "long-range" interaction are 
concerned, it was found that a nn repulsion U p d naturally 
disfavors the stripes while it does not affect the POLWC, 
where the charges are at larger distance. Nevertheless vi) 
a moderate U p d stabilizes the HFVS with respect to the 
FVS, thus strengthening the stability of this configura- 



tion. On the other hand, vii) the large-distance part of 
a truely long range interaction turns out to increase the 
energy of the POLWC thus favoring further the stability 
of an array of stripes. 

Finally, the effects of a coupling between holes and 
static lattice deformations were also investigated. In 
particular, the hopping was modulated by the breating 
displacement of oxygen ions in the j site Uj({uj}) — 
t p d ± auj, where the sign is positive if the Cu-0 distance 
is decreased, while it is negative if the Cu-0 distance in- 
creases. For simplicity, in this model we neglected the 
coupling between the 0-0 hopping and the lattice. In 
an analogous way the orbital energy on the copper site 
is modified as £i({ufc}) = £d + PJ2j(^ zU j)- Then the 
dimensionless hole-lattice coplings are A Q = a 2 /(Kt p d) 
and Xg — (3 2 / (Kt p d). At variance with the findings of 
Ref. nS, where stripe textures were only found as local 
minima, the slave-boson treatment produces absolutely 
stable stripe solutions. However, owing to the specific 
choice of lattice displacement, there occurs a frustration 
effect whenever holes (or electrons) are located on neigh- 
boring sites. This frustration of the lattice displacement 
is not present in the case of distant polarons. Therefore in 
stripe configurations (particularly in the filled ones), viii) 
the electron-phonon coupling turns out to be detrimental 
for stripes in comparison to the POLWC. Nevertheless, 
the absolute stability of HFVS at hole-doping x = 1/8 
survives in the presence of a moderate coupling. 

For concreteness, we report in Table I the energy dif- 
ference between the HFVS or the FDS and the POLWC 
configurations at x = 1/8 for a typical set of parameters 
in a 9 x 8 lattice. The HFVS stabilization effect of the 
long-range (LR) interaction is shown in the second line, 
while the destabilization due to the electron-phonon (EP) 
coupling is shown in the third line. The consideration of 
all the different interactions leads to a stable HFVS con- 
figuration as indicated in the fourth line. As shown in 
the last column, the FDS are instead always disfavored 
with respect to the POLWC. 
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Table I: Energy differences (per doped hole) between two 
different stripe configurations (HFVS and FDS) and the 
POLWC in a 9 x 8 (or 8 x 8) lattice with periodic bound- 
ary conditions. According to Ref. |io[ the parameter values 
are t pd = 1.3eV, t pp = 0.65eV, U d = 10.5eV, U p = 4eV. The 
LR coupling is such that U p d = 1.2eV with a screening length 
As = 32 lattice units. The EP couplings are X a — 0.5 and 

- Electron doping. While the above analysis was 
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mostly carried out for hole doping (and particularly for 
x = 1/8), it is quite important to extend the investiga- 
tion to electron-doped systems. While the general va- 
lidity of the above results is mantained, some important 
differences are also found. Specifically it was noticed that 
the charge in electron doped stripes is more localized. A 
comparison of the charge profile of hole-doped (HD) and 
electron-doped (ED) stripes is shown in Fig. 1, where the 
narrower charge distribution of ED stripes is apparent. 




FIG. 1. Charge density distribution in the (1,0) direction 
for FDS (i.e. in the (1,1) direction, above) and for HFVS (i.e. 
in the (0,1) direction, below) for both ED and HD systems. 
The parameters are the same as in Table I with x = ±1/8, 
but with A a = A/3 = 0. 

This effect is quite natural since the doped holes mostly 
reside on oxygen sites where due the smaller interaction 
and to the direct t pp hopping acquire a high mobility, 
while the doped electrons mostly reside on copper sites. 
Therefore the doped-electron mobility is lower because it 
only occurs via an intermediate charge transfer to oxygen 
sites. The smaller mobility and the consequent relative 
narrowness of the stripes in ED systems has two main 
consequences. The first is that the stripe formation is 
generically more difficult. Nevertheless stripe textures 
still remain the ground state for realistic values of the 
parameters. In particular, even in the presence of a siz- 
able electron-lattice coupling (which, as discussed above, 
favors the POLWC solution more than the stripes) sta- 
ble FDS textures are found. A second important effect is 
that FDS (with one doped electron per stripe site) become 
more stable than the HFVS. This latter finding suggests 



that ED cuprates should be more similar to the nick- 
elates, where diagonal stripes with one doped hole per 
site are detected. 

- Discussion. The above reported findings have im- 
portant consequences. First of all, the fact that FDS 
accomodate one doped charge per site, rather than one- 
half like the HFVS, increases the (average) distance be- 
tween the stripes. For instance, at |a;| = 1/8, in the FDS 
case the stripes are at distance oVds = 8a/y/2, while 
the HFVS are separated by dnFVS — 4a. As a con- 
sequence the interstripe antiferromagnetically correlated 
background is better preserved when FDS are formed. 
This effect is emphasized by the narrow width of the 
ED stripes, which is naturally accompanied by a lower 
transversal mobility of the stripes. It should also be no- 
ticed, that static disorder effects due to pinning of the 
stripes should comparatively be stronger in ED materi- 
als, where the stripes are expected to be less mobile. 

The formation of ordered arrays of stripes is nat- 
urally accompanied by an (incommensurate) spin re- 
sponse, which appears in neutron scattering at mo- 
menta displaced by an amount e from the AF wavevector 
Qaf = (±7r/a, ±7r/a). While this effect has indeed been 
observed in metallic HD materials, with s a±.x for x up 
to 1/8 and then saturating (e ~ constant :jE3, compara- 
tively detailed neutron scattering data are not available 
for ED systems. Our findings suggest that the incom- 
mensurate momenta should scale as £ ~ x/y/2. More- 
over, while in HD materials the spin excitation peaks ap- 
pear in the (1,0) and (0,1) directions consistently with 
the vertical ordering of the stripes, in ED materials the 
corresponding peaks should be observed in the (1,1) and 
(1, —1) directions, according to our finding of diagonal 
stripe formation. However, the substantial disorder ef- 
fects due to pinning could broaden the spin peaks making 
the observation of incommensurability more difficult. In 
this respect more local probesJike the NQR (successfully 
applied in the HD materialsO) can turn out to be very 
useful. 

As far as transport is concerned, no theory is presently 
available in the case of scattering due to (dynamically 
fluctuating) stripes. Therefore no definite statement can 
be put forward beside the following generic remarks. The 
first generic observation is that the narrower charge dis- 
tribution of stripes in ED systems, indicates a smaller 
transversal mobility. This, together with the larger sep- 
aration of the filled ED stripes accounts for the more ex- 
tended insulating region of the ED systems with respect 
to the HD ones. Specifically this could explain why the 
antiferromagnetic phase is quite more robust upon elec- 
tron doping than upon hole doping stressing the analogy 
between the stripes in ED systems and in nickelates. 

As far as the longitudinal (i.e. along the stripes) metal- 
licity is concerned, it can be noticed that HFVS should 
be more metallic than FDS. This is expected because 
in the longitudinal direction HFVS are similar to an ar- 
ray of quarter-filled onedimensional wires. If these wires 



3 



escape the longitudinal CDW commensurate instability, 
this ID structures should naturally be moreonetallic than 
ID chains with an empty or a full banc£3. If this is 
the case, a mechanisra-ifor superconductivity based on 
coupled metallic chainsa seems hardly applicable to the 
ED case at least at the low and intermediate doping 
regimes where (according to our calculations) the stripes 
are expected to be less mobile in the transverse direc- 
tion, nearly insulating in the longitudinal direction, and 
far one from the other. Of course this objection no longer 
holds when, by increasing doping, the metallic character 
of the system increases and our calculational scheme is 
no longer applicable. In this regime the stripes will be 
closer and likely far from the completely filled condition. 
In this latter doping region, however, the proximity to the 
quantum critical point observecLnear x = 0.1713 opens the 
way to alternative mechanismsB based on the scattering 
from critical charge fluctuations. 

The metallic properties of ED stripes should also be 
strongly affected by disorder effects, which, as already 
mentioned above, are expected to be more important 
when the stripes-are less mobile. In this regard, the re- 
cent observations that small amounts of interstitial api- 
cal oxygen suppress superconductivity (but not metallic- 
ity) giving rise to pseudogap effects is of obvious rele- 
vance. A natural speculation could be that the intersti- 
tial oxygen triggers a pinning of the stripes. The sup- 
pression of stripe fluctuations then induces the suppres- 
sion of superconductivity giving rise to a metallic (non- 
superconducting) state with disordered static (and lo- 
cally well-formed) stripes with the consequent stronger 
reduction of density of states near the Fermi level. 

Another generic consequence is that the diagonal 
stripes of the ED systems, should display rather different 
commensurability effects than in the HD materials. This 
is a consequence of the different spacing and ordering di- 
rection as well as of the different T' crystal structure of 
the ED cuprates. In particular, in the-ED case, the 1/8 
effects are not expected to be markedci since the stripes 
are quite far apart and the collective nature of the or- 
dered state is likely affected by disorder pinning. 

In conclusion, the occurrence of stripes in a rather re- 
alistic model for ED cuprates is a quite intriguing find- 
ing. According to our results, the specific features of 
the stripes in ED material are different from the fea- 
tures of stripes in HD ones. Therefore, if the occur- 
rence of stripes in electron-doped systems will be ex- 
perimentally confirmed, the comparison with the stripes 
already observed in the hole-doped materials will likely 
shed light on the physical mechanisms acting in both 
ED and HD cuprates. This calls for a search of experi- 
mental signatures and characterization of stripes in ED 
cuprates, hopefully leading to a better understanding of 
the cuprates. 
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